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Abstract 

The Focus constraint expresses the notion that so- 
lutions are concentrated. In practice, this constraint 
suffers from the rigidity of its semantics. To tackle 
this issue, we propose three generalizations of the 
Focus constraint. We provide for each one a com- 
plete filtering algorithm as well as discussing de- 
compositions. 



i-H 1 Introduction 



Many discrete optimization problems have constraints on the 
objective function. Being able to represent such constraints 
is fundamental to deal with many real world industrial prob- 
lems. Constraint programming is a promising approach to 
express and filter such constraints. In particular, several 
constraints have been proposed for obtaini ng well-balanced 
solutions iPesant and Regin, 2005; Scha us et al., 2007 
Petit and Regin, 20111. Recently, the Focus constraint [Pe- 
tit, 2012J was introduced to express the opposite notion. It 
captures the concept of concentrating the high values in a se- 
quence of variables to a small number of intervals. We recall 
its definition. Throughout this paper, X = [xq, a;i, . . . , Xn-i] 
is a sequence of variables and Sij is a sequence of indices of 



consecutive variables in X, such that s; 



\i,i + l,. 



>■ < i < j < n. We let \E\ be the size of a collection E. 

X 
J3 



Definition 1 (I Petit, 2012 1). Let yc be a variable. Let k and 
len be two integers, 1 < len < \X\. An instantiation of X U 
{yc} satisfies FocusfX, yc, len, k) iff there exists a set Sx 
of disjoint sequences of indices Si,j such that three conditions 
are all satisfied: (1) \Sx 
3s,;. ,- G Sx such that I € s 



'jj 



«j 



< yc (2) Wxi e X, XI > k ^ 
(3) Vsi J- G Sx, i — « + 1 < len 



Focus can be used in various contexts including cumu- 
lative scheduling problems where some e xcesses of ca pac- 
ity can be tolerated to obtain a solution [ Petit, 2012) . In 
a cumulative scheduling problem, we are scheduling activi- 
ties, and each activity consumes a certain amount of some re- 
source. The total quantity of the resource available is limited 
by a capacity. Exces ses can be represented by variables |De| 
Clercq et al., 20111. In practice, excesses might be toler- 



to its duration: it is cheaper to rent a machine during a sin- 
gle interval than to make several rentals. On the other hand, 
rental intervals have generally a maximum possible duration. 
Focus can be set to concentrate (non null) excesses in a 
small number of intervals, each of length at most len. 

Unfortunately, the usefulness of FOCUS is hindered by the 
rigidity of its semantics. For example, we might be able to 
rent a machine from Monday to Sunday but not use it on Fri- 
day. It is a pity to miss such a solution with a smaller number 
of rental intervals because Focus imposes that all the vari- 
ables within each rental interval take a high value. Moreover, 
a solution with one rental interval of two days is better than 
a solution with a rental interval of four days. Unfortunately, 
Focus only considers the number of disjoint sequences, and 
does not consider their length. 

We tackle those issues here by means of three generaliza- 
tions of Focus. SpringyFocus tolerates within each se- 
quence in Sij e Sx some values v < k. To keep the se- 
mantics of grouping high values, their number is limited in 
each Si J by an integer argument. WeightedFocus adds a 
variable to count the length of sequences, equal to the num- 
ber of variables taking a value v > k. The most generic 
one, WeightedSpringyFocus, combines the semantics 
of SpringyFocus and WeightedFocus. Propagation of 
constraints like these complementar y to an objective func- 
tion is well-know n to be important I Petit and Poder, 2008 1 
Schaus et al. , 2009 1 . We present and experiment with filtering 



ated by, for example, renting a new machine to produce more 
resource. Suppose the rental price decreases proportionally 



algorithms and decompositions therefore for each constraint. 

2 Springy FOCUS 

In DefinitionfT] each sequence in Sx contains exclusively val- 
ues V > k. In many practical cases, this property is too strong. 
Consider one simple instance of the problem in the introduc- 
tion, in Figure [T] where one variable Xi e X is defined per 
point in time i (e.g., one day), to represent excesses of capac- 
ity. Inintialy, 4 activities are fixed and one activity a remains 
to be scheduled (drawing A), of duration 5 and that can start 
from day 1 to day 5. If Focus(X, yc = 1, 5, 0) is imposed 
then a must start at day 1 (solution B). We have one 5 day 
rental interval. Assume now that the new machine may not 
be used every day. Solution (C) gives one rental of 3 days 
instead of 5. Furthermore, if len = 4 the problem will have 
no solution using Focus, while this latter solution still ex- 
ists in practice. This is paradoxical, as relaxing the condition 
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Lemma 1 (initialization). p{xo,v<:) ~ //'min(xo) < k, 
and n + 1 otherwise; ps{xo,v<) = n + 1; p(xo,w>) = 
1 ;/ inax(xo) > k and n + 1 otherwise; plen^Xg) — 1 if 
inax(a;o) > k and otherwise; card ixn) = 0. 

Lemma 2 {p{xi,v<)). If nim(xi) < k then p{xi,v<) = 
min(p(x/_i,u<),p(a;;_i,u>)), else p{xi,v<) == n + 1. 

Lemma 3 (ps{xi,v<)). Ifimn{xi)>k, ps{xi,v<)~n + 1. 
Otherwise, if plen{xi^i) G {0, ten — 1, len\ V card (xi_i ) 
— h then ps{xi,v<,) = n + 1, else ps{xi,v<,) = 
min(p5(a;i_i,u<),p(a;j_i,u>)). 

Lemma4 (p(x/, w>)). Ifin.sx.{xi) < k thenp{xi,Vy)='n+l. 
Otherwise, If plen{xi^i) G {0, Zen}, p{xi,Vy) = 
min(p(x;_i,u>) + l,p(a;/_i, u<) + 1), else p{xuVy) 
= min(p(a;;_i,w>),ps(a;;_i,t;<),p(x/_i,u<) + 1). 

Proposition 1 (plen{xi)). (by construction) If min 

{Ps{xi-i,v<),p{xi~i,v>))<p{xi-i,v<)+lA plen {xi_i)< 
ten then plen{xi) — plen{xi^i) + 1. Otherwise, if 
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Figure 1: (A) Problem with 4 fixed activities and one activity of 
lengtti 5 that can start from time 1 to 5. (B) Solution satisfying 
F0CUS(X, [1, 1], 5, 0), with a new machine rented for 5 days. (C) 
Practical solution violating Focus(X, [1, 1], 5, 0), with a new ma- 
chine rented for 3 days but not used on the second day. 

that sequences in the set Sx of Definition [T] take only values 
V > k deteriorates the concentration power of the constraint. 
Therefore, we propose a soft relaxation of Focus, where at 
most h values less than k are tolerated within each sequence 
inS'x- 

Definition 2. Let yc be a variable and k, len, h be three 
integers, 1 < len < |^|, < h < len — 1. An instantiation of 
X U {yc} satisfies SPRINGYFOCUSfX, y^ len, h, k) iff there 
exists a set Sx of disjoint sequences of indices Si,j such that 
four conditions are all satisfied: (1) l^xl < yc (2) Vx; G X, p{xi,Vy)) < n + 1 then plen{xi) = 1, else plen{xi) = 0. 

Proposition 2 (card (x/)). (by construction) If plen{xi) — 1 
then card jxi ) — 0. Otherwise, if p{xi , w> ) = n + 1 then 

card (xi) — cardixi^i) + 1, else card (xi) — card (xi^-\ ). 
Bounds consistency (BC) on SpringyFocus is equiva- 
lent to domain consistency: any solution can be turned into 
a solution that only uses the lower bound min(a;/) or the up- 
per bound max(a;;) of the domain Dj xi) of each xi G X 
(this observation was made for Focus I Petit, 2012] ). Thus, 
we propose a BC algorithm. The first step is to traverse X 
from xq to 2:„_i, to compute the minimum possible number 
of disjoint sequences in Sx (a lower bound for yc), the fo- 
cus cardinality, denoted fc{X). We use the same notation for 
subsequences of X. fc{X) depends on fc, len and h. 

Definition 3. Given xi G X, we consider three quantities. 
(1) p{xi,v<,) is the focus cardinality of [xo,xi, ... ,xi], as- 
suming xi < k, and Vsi^j G S[^g^^^^,„^^,],j ^ I. (2) 
Ps{xi,v<:) is the focus cardinality of [xq, Xi, . . . ,Xi], assum- 
ing xi < kand3i,St^i G ^'[j.,,^^,^^...^^,,]. (3) p{xi,Vy) is the 
focus cardinality of [xo,xi, ... , xi] assuming xi > k. 

Any quantity is equal to n -\- 1 if the domain D{xi) of xi 
makes not possible the considered assumption. 

Property 1. ps{xo,v<) = ps(x„_i,i;<) = n + 1, and 
fc(X) = min(p(x„_i,t;<),p(x„_i,i;>)). 



Algorithm 1: MinCards(X, len, k, hy. integer mamx 

1 pre :— new Integer[|X|] [4] [] ; 

2 for / G 0..n — 1 do 

3 pre[^][0] :— new Integer[2]; 

4 for j e 1..3do pre[/][j] :— new Integer [1]; 

5 Initialization Lemmani; 

6 for / tE l..n — 1 doT^emmaspj Isl pland Propositions [lland pj; 

7 return pre; 



Proof By construction from Definitions l2] and l3] 
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To compute the quantities of Definitionl3]for xiEX we use 
plen{xi), the minimum length of a sequence in S[xo.xi,...,xi] 
containing xi among instantiations of [xq, xi, . . . ,xi] where 
the number of sequences is fc{[xQ,xi, . . . , xi]). plen{xi)=Q 

ifVsj j G S[rcg^xi....,xi]yj ¥" I- card (xi) is the minimum num- 
ber of values v < k in the current sequence in -S'fj.Q ^.^ j;,], 
equal to if Vsi.j G •Sj^j^j.^ x^],j ^ I. card (xi) assumes 
that Xl > k. It has to be decreased it by one if xi < k. For 
sake of space, proofs of next lemmas are given in Appendix. 



Algorithm [T] implements the lemmas with pre [I] [0] [0] = 
p{xi,v<), pre[l][0][l] = ps^(xi,v<), pre[l][l] = p{xi,Vy), 
pre[l][2] = plen{xi), pre[l][3] = card (xi). 

The principle of Algorithmic] is the following. First, lb ~ 
fc{X) is computed with x„_i. We execute Algorithmfllfrom 
Xq to a;„_i and conversely (arrays pre and suf). We thus 
have for each quantity two values for each variable xi. To ag- 
gregate them, we implement regret mechanisms directly de- 
rived from Propositions I2] and [T] according to the parameters 
len and h. Line 4 is optional but it avoids some work 

when the variab le yc is fix ed, thanks to the same property 
as Focus (see [Petit, 20121). Algorithm [2] performs a con- 
stant number of traversals of the set X. Its time complexity 
is 0{n), which is optimal. 

3 Weighted FOCUS 

We present WeightedFocus, that extends Focus with a 
variable Zc limiting the the sum of lengths of all the sequences 
in Sx, i.e., the number of variables covered by a sequence in 
Sx. It distinguishes between solutions that are equivalent 
with respect to the number of sequences in Sx but not with 
respect to their length, as Figure|2]shows. 



Algorithm 2: FilteriNG(X, y^, len, k, h): Set of variables 

1 pre := M[NCards(X, len, fc, ?i) ; 

2 Integer Ih -.— min(pre[n — 1][0] [0], pre[n — 1][1]); 

3 if min(ac) < ibthen -D(i/c) := -D(yc) \ [min(i/c), ifc[; 

4 if min(yc) — inax(yc) then 

5 su/ :— MinCards([xti-i, Xn-2, ■ ■ ■ -Z^^]; ^en, /c, /i) ; 

6 for / e 0..n — 1 do 

7 if pre[;][0][0] + su/[n - 1 - i][0][0] > max(i/<;) tlien 

8 Integer regret :— 0; Integer add :— 0; 

9 if pre[i][l] < pre[;][0][l] tJien add := add + 1; 
10 if suf[n - 1 - /][1] < suf[n - 1 - i][0][l] tlien 

add:^add+l; 
u if pre[;][2] + su/[n - 1 - l][2] - 1 < len A 

pr-e[;][3] + suf[n - 1 - /][3] + add - 1 < /i then regret 
:= 1; 

12 if 

pre[;][0][l] + suf[n - 1 - i][0][l] - regret > inax(j/c) 
then £1(2;;) := -D(a;i)\ [min(a:i), fc]; 

13 Integer regret :— 0; 

14 if pre[i][2] + s«/[n - 1 - l][2] - 1 < kn A 
pre[Z][3] + SM/[n - 1 - i] [3] - 1 < h then regret := 1; 

15 if pre[^][l] + 5w/[n — 1 — /][!] — regret > inax(j/c) then 

16 D{xi) :— D{xi)\]k,uiR-i^{xi)]\ 

17 return X U { yc } ; 



Definition 4. Le/ j/c and Zc be two integer variables 
and k, len be two integers, such that 1 < len < 
\X\. An instantiation of X U {i/c} U {^c} satisfies 
WEIGHTEDFocusfX, j/c, len, k, Zc) iff there exists a set Sx 
of disjoint sequences of indices Si,j such that four conditions 
are all satisfied: (1) \Sx\ < Vc (2) Vx; ^ X, xi > k ^ 
3si j G Sx such that I € Si j (3) Vs^ j € Sx, i — * + 1 < ^^^ 

Definition 5 (| |Petit, 2012| ). Given an integer k, a variable 
xi ^ X is: Penalizing, (Pk), iff rninixi) > k. Neutral, 
(Nk), iff max{xi) < k. Undetermined, (Uk), otherwise. We 
say xi G Pk iffxi is labeled Pk, and similarly for Uk and Nk- 
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Figure 2: (A) Problem with 4 fixed activities and one activity of 
lengtii 5 that can start from time 3 to 5. We assume D{yc) = {2}, 
len = 3 and fc = 0. (B) Solution satisfying WeightedFocus with 
Zc = 4. (C) Solution satisfying WeightedFocus with Zc = 2. 



Dynamic Programming (DP) Principle Given a partial in- 
stantiation Ix of X and a set of sequences Sx that cov- 
ers all penalizing variables in Ix, we consider two terms: 
the number of variables in Pk and the number of undeter- 
mined variables, in Uk, covered by Sx- We want to find a 
set Sx that minimizes the second term. Given a sequence 
of variables Sij, the cost cst{si,j) is defined as cst{sij) — 
{p\xp £ Uk,Xp G Sij}. We denote cost of Sx, cst{Sx), 
the sum cst{Sx) = J2si eSx cst{sij). Given Ix we con- 
sider |Pfe| = \{x, G K}|.'' We have: Es,.,es\^u\ = 



We start with explaining the main difficulty in building a 
propagator for WeightedFocus. The constraint has two 
optimization variables in its scope and we might not have a 
solution that optimizes both variables simultaneously. 

Example 1. Consider the set X — [a;o, 2:1, . . . ,^5] 
with domains [1, {0, 1}, 1, 1, {0, 1}, 1] and 

WeightedFocus (X, [2, 3], 3, 0, [0, 6]), solution 

Sx = {so,2,S3,5}> Zc = 6, minimizes i/c = 2, while 
solution Sx = {so,i, ^2,3, 55,5}. Vc = 3, minimizes Zc = 4. 

Examplefllsuggests that we need to fix one of the two opti- 
mization variables and only optimize t he other one. Our algo - 
rithm is based on a dynamic program | Dasgupta et ah, 2006) . 
For each prefix of variables [xq, xi, . . . , a;jj and given a cost 
value c, it computes a cover of focus cardinality, denoted 
Sc.j, which covers all penalized variables in [xq,xi, . . . ,Xj\ 
and has cost exactly c. If Sc.j does not exist we assume that 
Sc.j = 00. Sc.j is not unique as Examplel2]demonstrates. 

Example 2. Consider X = [xq,xi,. . . ,xt\ and 
WeightedFocus(A:, [2, 2], 5, 0, [7, 7]), with D{xi) = {!}, 
i e 1,1 = {0,2,3,5,7} and D{xi) = {0,1}, 
i G {0,1,... 7} \ /. Consider the subsequence of 
variables [xo, . . . ,2:5] and S'l^s. There are several sets of 
minimum cardinality that cover all penalized variables in the 
prefix [a;o, ■ • • , a^s] and has cost 2, e.g. S\^ = {so,2, S3, 5} 
or Sfj^ — {50,4,55,5}- Assume we sort sequences by their 
starting points in each set. We note that the second set is 
better if we want to extend the last sequence in this set as the 
length of the last sequence 55^5 is shorter compared to the 
length of the last sequence in S\ g, which is 53^5. 

Example [2] suggests that we need to put additional condi- 
tions on Sc^j to take into account that some sets are better than 
others. We can safely assume that none of the sequences in 
Sc,j starts at undetermined variables as we can always set it 
to zero. Hence, we introduce a notion of an ordering between 
sets Sc,j and define conditions that this set has to satisfy. 

Ordering of sequences in Sc.j. We introduce an order 
over sequences in Sc.j ■ Given a set of sequences in Scj we 
sort them by their starting points. We denote last{Sc,j) the 
last sequence in Sc.j in this order If Xj G last{Sc.j) then 
\last{Scj)\ is, naturally, the length of last{Sc.j), otherwise 
\last{Sc,j)\ = 00. 

Ordering of sets Sc.j, c G [0, max(zc)], j G {0, 1, . . . , n — 
1}. We define a comparison operation between two sets S, 



and Sc'. 



< Sc 



iff 15-^ 



c,3 



< ISc'^A or \Sc.,\ = 
J,). Note that we do not 



j' "i \>^C,J 

\Sc'j'\ Sind last{Sc,j) < last{S, 
take account of cost in the comparison as the current defini- 
tion is sufficient for us. Using this operation, we can compare 
all sets Sc,j and S'c of the same cost for a prefix [xq, . . . , Xj]. 
We say that Scj is optimal iff satisfies the following 4 condi- 
tions. 

Proposition 3 (Conditions on S^j). 



Sc,j covers all Pk variables in 



a;o,a;i. 



CSt{Sc,j) = c, 
Vs,i_g G Sc,J,Xh 



Uk, 



Sc.j is the first set in the order among all sets that satisfy 
conditionsU]^ 



Xj-1 Xj in Pk Xj-i Xj in Uk Xj-i x, in Nk 
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Figure 3: Representation of one step of Algorithm |3] 

As can be seen from definitions above, given a subse- 

Sc,j is not unique and might 
?c'j'|' c = c' and j == /, 



quence of variables xa, . . . . 

not exist. However, if \Sc.j\ — \S, 

then last{Sc.j) — last{Sc'.j')- 

Example 3. Consider WeightedFocus from Example |2] 
Consider the subsequence [a;o,a;i]. 5o.i — {so,o}. "S*!,! = 
{so,i}- Note that 6*2,1 does not exist. Consider the sub- 
sequence [xq, ... ,1^]. We have Sa.r^ = {so,o, S2,3, 55,5}, 
-Si, 5 = {so,4,S5,5} a«^ 'S'2,5 = {50,3,55,5}- By definition, 
last{So^5) = 55,5, last^Si^^) = 55,5 and last{S2,5) = s^^^. 
Consider the set Si^^. Note that there exists another set 
^'1 5 — {sofi, 82,5} that satisfies conditions\lH3\ Hence, it 
has the same cardinality as Si 5 and the same cost. However, 
S'1,5 < Si.5 as \last{Si^5)\ = 1 < \last{S[^5)\ = 3. 

Bounds disentailment Each cell in the dynamic program- 



ming table fcj, c e [0, z^], j G {0, 1 



,u 



1}, where 



\Pk\, is a pair of values Qcj and Icj, 



z^ = max{zc) 

fc,j — {ic.jjc.j}, stores information about Scj. Namely, 
1c,j — \Sc,j\, lc,j — \lo,st{Sc,j)\ if last{Sc,j) 7^ 00 and 00 
otherwise. We say that fcj/qcj/lcj is a dummy (takes a 
dummy value) iff /c,j — {oo,oo}/gc.j = ^x)/lcj — 00. If 
2/1 = 00 and y2 — 00 then we assume that they are equal. We 
introduce a dummy variable a:_i, D{x^i) = {0} and a row 
f-i.j, j = — 1, • • . , n — 1 to keep uniform notations. 



Algorithm 3: Weighted Focus(a;o, 



5 -^n—l} 



1 for c G — 1 . . 2 j^ do 

2 for J G — l..n — 1 do 

3 /c,j •*- {00,00}; 

4 /o.„i <- {0,0}; 

5 for J G 0..n — 1 do 

6 fore g 0..J do 

7 \ixj g P^ then /* penalizing */ 

8 if (/c.j-i £ [1, len)) V (gc.j-i — oo) then 
« /cj ^ {qc,j-i, la. j-i + 1}; 

10 else /cj *- {9c. 3-1 + 1, 1}; 

11 \ixj g C/fc then /* undetermined */ 

12 if(;c-l,j-l e [l,ien) A qc-l,j-l = 
(?c,j-i) V (gc,j-i — 00) then 

/cj <- {<3c-i,3-i, 'c-i,j-i + 1} else 
/c.j •*- {?c,3-i,oo} 

13 if a:;j G A'';^ then /* neutral */ 
IJ /cj <- {9c, 3 -1,00} 



IS return /; 



Algorithm [3] gives pseudocode for the propagator. The in- 
tuition behind the algorithm is as follows. We emphasize 
again that by cost we mean the number of covered variables 
inC/fe. 

If Xj e Pk then we do not increase the cost of 5c j com- 
pared to Sc,j^i as the cost only depends on Xj € Uk. Hence, 
the best move for us is to extend last{Sc.j-i) or start a new 
sequence if it is possible. This is encoded in hnes|9]and 10 of 
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D(x^) 
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D(X2) 
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{1,00} 
{1.5} 


{3,1} 
{2,1} 
{2,1} 


{3,00} 
{2,00} 
{2,2} 


{4,1} 
{3,1} 
{2,3} 



Table 1: An execution of Algorithm[3]on WeightedFocus 
from Examplel2] Dummy values fcj are removed. 

the algorithm. Figure|3ja) gives a schematic representation of 
these arguments. 

If Xj e Uk then we have two options. We can obtain 
Sc.j from S'c-i.j-i by increasing cst{Sc-i.j-i) by one. This 
means that Xi will be covered by last{Scj). Alternatively, 
from Scj-i by interrupting last{Sc.j^i). This is encoded in 
line[T2]of the algorithm (Figure[3]^b)). 

If Xj E Nk then we do not increase the cost of Scj com- 
pared to Sc.j-i. Moreover, we must interrupt last{Sc,j-i), 
line 14 (Figurel3jc), ignore the gray arc). 

First we prove a property of the dynamic programming 
table. We define a comparison operation between fcj and 
fc'.j' induced by a comparison operation between Scj and 
Sc',]'- Ic.j < fc',]' if {qc,] < qc'.j') or (qcj = q^ ,y and 
Ic.j < Ic'.j')- In other words, as in a comparison operation 



between sets, we compare by the cardinality of sequences, 
I Scj I and I Sc' ,j'\, and, then by the length of the last sequence 
in each set, ^asi(S'c.j) and last{Sc\j'). We omit proofs of the 
next two lemmas due to space limitations (see Appendix). 

Lemma 5. Consider WeightedFoc\JS{X, yc, len, k, z^). 
Let f be dynamic programming table returned by Algo- 
rithm [5] Non-dummy elements fcj are monotonically non- 
increasmg in each column, so that fc'j < fcj, < c < c' < 
z^,j^[0,...,n^l]. 

Lemma 6. Consider WeightedFocus (X, j/c, ^era, A;, Zc). 
The dynamic programming table fcj — {icj^^cj} c € 



[0, z,"], J = 0, . . 



1, is correct in the sense that if fcj 



exists and it is non-dummy then a corresponding set of se- 
quences Sc.j exists and satisfies conditions |7}p] The time 
complexity of Algorithmpjis 0{ninax{zc)). 

Example 4. Table[I\shows an execution of Algorithm^on 
WeightedFocus /rom iixampte [2] Note that \Po\ = 5. 



Hence, 



,u _ 



max{zc) — I -Pol = 2. As can be seen from 
the table, the constraint has a solution as there exists a set 
82,7 = {50,3, ss,?} such that \S2j\ = 2. 



Bounds consistency To enforce BC on variables x, we 
compute an additional DP table h, bcj, c G [0, z^], j G 
[— 1, n — 1] on the reverse sequence of variables x. 

Lemma 7. Consider WeightedFocus(X, j/c, /en, fc, Zc)- 
Bounds consistency can be enforced in 0{n niax(zc)) time. 

Proof. (Sketch) We build dynamic programming tables / and 
b. We will show that to check if Xi — v has a support it 
is sufficient to examine 0{z^) pairs of values fci,i-i and 
6c2 ,n-j-2j Ci , C2 G [0, z^] which are neighbor columns to the 
ith column. It is easy to show that if we consider all possible 
pairs of elements in fci,i-i and 6c2.n-i-2 then we determine 



if there exists a support for Xi = v. There are 0(2^ 



') 



such pairs. The main part of the proof shows that it sufficient 
to consider 0{z^) such pairs. In particular, to check a sup- 
port for a variable-value pair x, ~ v, v > k, for each f^^i-i 
it is sufficient to consider only one element 6c2.n-i-2 such 



that br 



_2 is non-dummy and C2 is the maximum value 



that satisfies inequality ci + C2 + 1 < z^ . To check a support 
for a variable-value pair Xi = v, v < k, for each fci.i-i it 
is sufficient to consider only one element 6c2,n-J-2 such that 
bc2.n-i-2 is non-dummy and C2 is the maximum value that 
satisfies inequality ci + C2 < z^ . D 

We observe a useful property of the constraint. If there 
exists /c,n-i such that c < max{zc) and qc,n-i < rnax{yc) 
then the constraint is BC. This follows from the observation 
that given a solution of the constraint Sx, changing a variable 
value can increase cst{Sx) and \Sx\ by at most one. 

Alternatively we can decompose WeightedFocus using 
0{n) additional variables and constraints. 

Proposition 4. Given FocusfX, yc, len, k), let Zc be a vari- 
able and _B=[6o, 61, ■ ■ • , &n-i] be a set of variables such that 
VbieB, _D(6()={0, 1}. WEIGHTEDPocusfX, j/c, len, k, Zc) 
^ FocusfX, yc, len, k) A \\/l, < I < n, [{xi < k) A [h = 
0)] V [{xi >k)A {hi - 1)]] A Eie{o,i,...,n-i} bi < z,,. 

Enforcing BC on each constraint of the decomposition is 
weaker than BC on WeightedFocus. Given i; e X, a 
value may have a unique support for Focus which violates 
Sielo 1 n-i> ^' — ■^c, and conversely. Consider n=h, 
Xq=X2—'^, 2:3=0, and D{x\)—D(x^—{^,\'\, yc—^, Zc—3, 
k—0 and len=3. Value 1 for X4 corresponds to this case. 



4 Weighted Springy FOCUS 

We consider a further generalization of the FOCUS constraint 
that combines SpringyFocus and WEIGHTEDFOCUS. We 
prove that we can propagate this constraint in 0{n max(zc)) 
time, which is same as enforcing BC on WeightedFocus. 

Definition 6. Let yc and Zc be two variables and k, len, 
h be three integers, such that 1 < len < \X\ and < 
h < len — 1. An instantiation of X U {yc} U Zc satisfies 
WElGHTEDSPRlNGYFocusfX, yc, len, h, k, Zc) iff there ex- 
ists a set Sx of disjoint sequences of indices Sij such that 
five conditions are all satisfied: (1) \Sx\ < yc (Z) Vx/ € X, 
Xi > k => 3si,j G Sx such that I e Sij (3) Vs^j- £ Sx, 
\{l G Sij, xi < k}\ < h (4) Vsjj G Sx, j - i + 1 < len, 



> k and Xj > k. (5) ^^ 



i&Sj, 



^«J 



< Zc 



We can again partition cost of S into two terms. 
Es,,,esl'5«jl = Es.,,GsCs*(si,i) + l^fel- However, 
cst{si,j) is the number of undetermined and neutral variables 
covered Sij, cst{sij) — {p\xp £ Uk^ Nk,Xp G Si^-} as we 
allow to cover up to h neutral variables. 

The propagator is again based on a dynamic program that 
for each prefix of variables [xq, xi, . . . , Xj] and given cost 
c computes a cover Sc.j of minimum cardinality that covers 
all penalized variables in the prefix [xo, xi, . . . , Xj\ and has 
cost exactly c. We face the same problem of how to com- 
pare two sets Si and S^ of minimum cardinality. The is- 
sue here is how to compare last{Slj) and last{Slj) if they 



cover a different number of neutral variables. Luckily, we 
can avoid this problem due to the following monotonicity 
property. If last{Sl A and last{Sl A are not equal to in- 
finity then they both end at the same position j. Hence, if 
last{Sc j) < last{Sc j) then the number of neutral variables 
covered by last{Sl ) is no larger than the number of neutral 
variables covered by last{S'^ ). Therefore, we can define or- 
der on sets Sc,j as we did in SectionlSlfor WeightedFocus. 
Our bounds disentailment detection algorithm for 
WeightedSpringyFocus mimics Algorithm [3] We omit 
the pseudocode due to space limitations but highlight two 
not-trivial differences between this algorithm and Algo- 
rithm I3] The first difference is that each cell in the dynamic 
programming table fc,j, c G [0, z^], j G {0, 1, . . . , n — 1}, 



where 



max{zc) 



^fe|, is a triple of values Qcj, 
{Qc,j,lc,j,h'c.j}- The new parameter 
he J stores the number of neutral variables covered by 
last{Sc.j)- The second difference is in the way we deal with 

iVfc then we have two options 
ij_i by increasing 



Ic.j and he J, jc.j 



neutral variables. If x, G 



now. 



We can obtain Sc,j from S, 



cst{Sc-i.j-i) by one and increasing the number of covered 
neutral variables by last{Sc,j~i) (Figure Islc), the gray arc). 
Alternatively, we can obtain Sc,j from Sc,j-i by interrupting 
last{Sc.j-i) (Figure Isle), the iDlack arc). BC can enforced 
using two modifications of the corresponding algorithm for 
WeightedFocus (a proof is given in Appendix). 

Lemma 8. Consider WeightedSpringyFocus(X, j/c, 
len, h, k, Zc). BC can be enforced in 0{n niax(zc)) time. 

WeightedSpringyFocus can be encoded using the 
cost-REGULAR constraint. The automaton needse 3 coun- 
ters to compute len, yc and h. Hence, the time complexity of 
this encoding is O(n^). This automaton is non-deterministic 
as on seeing v < k, it either covers the variable or interrupts 
the last sequence. Unfortunately the non-deterministic cost- 
REGULAR is not implemented in any constraint solver to our 
knowledge. In contrast, our algorithm takes just 0{n?) time. 
WeightedSpringyFocus can also be decomposed using 
the Gcc constraint iRegin, 19961. We define the following 
variables for alH G \{),max[yc) — l\ andj G [0,7i— 1]: Si the 
start of the ith sub-sequence. D{Si) = {Q, ..,n + max{yc)}', 
Ei the end of the ith sub-sequence. D{Ei) — {0, ..,n + 
max{yc)}', Tj the index of the subsequence in Sx contain- 
ing Xj. D{Tj) = {0, .., ■max{yc)}', Zj the index of the sub- 
sequence in Sx containing Xj s.t. the value of Xj is less 
than or equal to k. D{Zj) — {0, ..,max{yc)}', lastc the 
cardinality of Sx. D{lastc) ~ {0, ..,max{yc)}', Card, a 
vector of max{yc) variables having {0,..,^} as domains. 
WeightedSpringyFocus(X, yc, len, h, k, Zc) <^ 

{xj <k)y Zj=Q; {xj <k)y Tj>0; 

{xj > k)\/ {Tj = Zj); {Tj< lastc); 

{Tj / i) V U > S,_i); {T, / i) V {j < i5.-i); 

(j > lastc) V (T, = i)\/(j < S^-i) V {j > E,^i); 

Vg G [l,max{yc) — 1], g > lastc V Sq > Eq-i; 

Vg G [0, max{yc) - 1] g > lastc V B, > 5,; 

Vg G [0,max{yc) — 1] g > lastc V len > (Eg — Sq); 

lastc < Vc-, GccilTo, .., r„_i], {0}, [n - Zc]); 

Gcc([Zo, .., Zn-i],{l, ..,niax{yc)}. Card); 



Table 2; SLS with WeightedFocus and its decomposition. 
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5 Experiments 

We used the Choco-2. 1 .5 solver on an IntelXeon 2.27GHz for 
the first benchmarks and IntelXeon 3.20GHz for last ones, 
both under Linux. We compared the propagators (denoted by 
F) of WeightedFocus and WeightedSpringyFocus 
against two decompositions (denoted by Di and D2), using 
the same search strategies, on three different benchmarks. 
The first decomposition, restricted to WEIGHTEDFOCUS, is 
shown in propositionH] while the second one is shown in Sec- 
tion |4] In the tables, we report for each set the total number 
of solved instances (#n), then we average both the number of 
backtracks (#b) and the resolution time (T) in seconds. 
□ Sports league scheduling (SLS). We extend a single round- 
robin problem with n = 2p teams. Each week each team 
plays a game either at home or away. Each team plays ex- 
actly once all the other teams during a season. We minimize 
the number of breaks (a break for one team is two consecutive 
home or two consecutive away games), while fixed weights in 
{0, 1} are assigned to all games: games with weight 1 are im- 
portant for TV channels. The goal is to group consecutive 
weeks where at least one game is important (sum of weights 
> 0), to increase the price of TV broadcast packages. Pack- 
ages are limited to 5 weeks and should be as short as possible. 
Table |2] shows results with 16 and 20 teams, on sets of 50 in- 
stances with 10 random important games and a limit of 400K 
backtracks. max(?/c) = 3 and we search for one solution 
with h < 7 (instances n-1), h < 6 (n-2) and h < 5 (n-3). 
In our model, inverse-channeling and AllDifferent con- 
straints with the strongest propagation level express that each 
team plays once against each other team. We assign first the 
sum of breaks by team, then the breaks and places using the 
DomOverWDeg strategy. 

D Cumulative Scheduling with Rentals. Given a horizon of 
n days and a set of time intervals [si, Ci], i £ {1,2,... ,p}, 
a company needs to rent a machine between k and Ui times 
within each time interval [si,ei]. We assume that the cost 
of the rental period is proportional to its length. On top 
of this, each time the machine is rented we pay a fixed 
cost. The problem is then defined as a conjunction of one 
WeightedSpringyFocus(X, ?/c, ^en, /i,0,Zc) with a set 
of Among constraints. The goal is to build a schedule for 
rentals that satisfies all demand constraints and minimizes 
simultaneously the number of rental periods and their total 
length. We build a Pareto frontier over two cost variables, as 
Figure [4] shows for one of the instances of this problem. We 
generated instances having a fixed length of sub-sequences of 
size 20 (i.e., len = 20), 50% as a probability of posting an 
Among constraint for each {i,j) s.t. j > z + 5 in the sequence. 




Figure 4: Pareto frontier for Scheduling with Rentals. 



Table 3: Scheduling with Rentals. 
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Each set of instances corresponds to a unique sequence size 
({40, 43, 45, 47, 50}) and 20 different seeds. We summarize 
these tests in table l3] Results with decomposition are very 
poor. We therefore consider only the propagator in this case. 

□ Sorting Chords. We need to sort n distinct chords. Each 
chord is a set of at most p notes played simultaneously. The 
goal is to find an ordering that minimizes the number of notes 
changing between two con secutive cho rds. The full descrip- 
tion and a CP model is in |Petit, 201 21. The main difference 
here is that we build a Pareto frontier over two cost variables. 
We generated 4 sets of instances distinguished by the num- 
bers of chords ({14, 16, 18, 20}). We fixed the length of the 
subsequences and the maximum notes for all the sets then 
change the seed for each instance. 

Tables|2)[3]and|4]show that best results were obtained with 
our propagators (number of solved instances, average back- 
tracks and CPU time over all the solved instances^. FigureH 
confirms the gain of flexibility illustrated by Figure [Tlin Sec- 
tion l2j allowing h = 1 variable with a low cost value into 
each sequence leads to new solutions, with significantly lower 
values for the target variable j/c- 



6 Conclusion 

We have presented flexible tools for capturing the concept of 
concentrating costs. Our contribution highlights the expres- 
sive power of constraint programming, in comparison with 
other paradigms where such a concept would be very difficult 
to represent. Our experiments have demonstrated the effec- 
tiveness of the proposed new filtering algorithms. 



'while the technique that solves the largest number of instances 
(and thus some harder ones) should be penalized. 
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Appendix 



Proofs of Lemmas [T] to m ommit the obvious cases where 
quantities take the default value ?i + 1. 

A Proof of Lemma [1] 

From item 4 of Definition |2] a sequence in Sx cannot start 
with a value v < k. Thus, pg{xQ, t;<) — n + 1 and card ixn) 
— 0. If xq can take a value v > k then by Definition l3] 

p{xq, u>) = 1 and plen{xo) = 1. D 

B Proof of Lemma |2] 

If min(x/) < kthenps{xi^i,v<:) must not be considered: it 
would imply that a sequence in Sx ends by a value v < k for 
xi-^i. From Property [T] the focus cardinality of the previous 
sequence is min(p(x;_i, w<),p(a;;_i,t'>)). D 



k we have three cases to consider (1) If ei- 

i) = or plen{xi-i) — len then from item 



C Proof of Lemma |3] 

If u\va{xi) < 
ther plen(xi^ 

3 of Definition l2] a sequence in Sx cannot start with a value 
Vi < k: ps{xi,v<c) — n + 1. (2) If plen{xi^i) = len — 1 
then from Defiinition|2]the current variable Xi cannot end the 
sequence with a value Vi < k. (3) Otherwise, from item 3 of 
Definitionl2] p{xi-i,v<:) is not considered. Thus, from Prop- 
erty [I|ps(xi,t)< ) ^ mm{ps_{xi-i,v<),p{xi-i,v::,)). D 

D Proof of Lemma |4] 

If plen{xi^i) e {0, len} a new sequence has to be con- 
sidered: Ps{xi~itV<) must not be considered, from item 
3 of Definition l2| Thus, p{xi,vy) — min(p(a;i_i,?;>) + 
l,p(a;;_i, v<) + 1). Otherwise, either a new sequence has to 
be considered {p{xi^i , v< ) + 1) or the value is equal to the fo- 
cus cardinality of the previous sequence ending in x;_i. D 

E Proof of Lemma |5] 

First, we prove two technical results. 

Lemma 9. Consider WElGHTEDFocus([a;o, ■ • • ,2:n-i], 
yc, len, k, Zc). Let f be dynamic programming table returned 
by Algorithm[3\ Then the non-dummy values of fc.-i are con- 



secutive in each column 



so that there do not exist c, c . 
< c < c' < c" < 7JJ , such that /c 
/cji /c".j fl'"^ non-dummy. 



is dummy and 



Proof. We prove by induction on the length of the sequence. 
The base case is trivial as /o.-i = {0,0} and fc.-i = 
{oo,cx)}, c G [—1] U [1,2^^]. Suppose the statement holds 
for j — 1 variables. 

Suppose there exist c, c', c", < c < c' < c" < z^ , such 
that fc'.j is dummy and fcj, fc" ,j are non-dummy. 

Case 1. Consider the case Xj e Pk- By Algo- 
ithmls] lines [9] and 



h 



10 



9c,i G [gc,i-i,9cj-i + 1], qc\j G 



j-i,qc',j-i 

fc',j is dummy and /cj, /, 



1] andqc"j G [gc"j-i, gc"j-i + !]■ As 

c'j-i must 



are non-dummy, /, 



be dummy and fc,j-i, fc".j-i must be non-dummy. This vi- 
olates induction hypothesis. 

Case 2. Consider the case Xj £ Uk- By Algo- 
rithm Is] lines 12 and 12 Qcj — min{qc-ij-i,qcj-i), 
Qc'j = min{qc'-ij-i,qc',j-i) and g^'-j = 

min{qc"-i,j-i,qc".j-i)- As fc\j is dummy, then both 
fc'-i.j-i and fc'.j-i must be dummy. As fc.j is non-dummy, 
then one of /c-ij-i and fc.j-i is non-dummy. As /c",j 
is non-dummy, then one of fc"-i,j-i and fc".j-i is non- 
dummy. We know that c-1 < c < c' -1 < c' < c" -1 < c" 
or c < c' < c". This leads to violation of induction hypothe- 
sis. 

Case 3. Consider the case Xj E Nk- By Algorithm [3] 



line [T4j qcj = qcj-i, Qc'j = Qc/j-i and 9c"j — 
qc",j~i- Hence, /c'j-i is dummy and fc,j-i, fc",j-i are 
non-dummy. This leads to violation of induction hypothe- 
sis, n 

Propositions. Co««Wer WElGHTEDFocus([a:o, • ■ • ,a;„_i], 

yc, len, fc, Zc). Let f be dynamic programming table returned 
by Algorithm^ The elements of the first row are non-dummy: 
fo.j> j = ~^T ■ ■ ■ 1^ ore non-dummy. 

Proof. We prove by induction on the length of the sequence. 
The base case is trivial as /o.-i — {0, 0}. Suppose the state- 
ment holds for j — 1 variables. 

Case 1. Consider the case Xj e Pk. As /o.j-i is non- 
dummy then by AlgorithmlS] linesl9- 10 foj is non-dummy. 

Case 2. Consider the case Xj e (Tfc. Consider the condition 
(Li , 1 £ [IJen) A q-i.j-i = 9o,j-i) V (goj-i = oo) at 
line|l2l By the induction hypothesis, qoj-i ^ oo. By the 

oo. 
is 



i.ti 



initialization procedure of the dummy row, q- 
Hence, this condition does not hold and, by line 12 
non-dummy. 



/( 



Oj 



Case 3. Consider the case Xj G A^^. As /oj-i is non- 

is non-dummy. 

n 



dummy then by Algorithmic] lines 14 /o 



We can now prove Lemma l5] 

Proof By transitivity and consecutivity of non-dummy val- 
ues (Lemma[9]) and the result that all elements in the 0th row 
are non-dummy (Proposition|5]l, it is sufficient to consider the 
case c' — c+1. 

We prove by induction on the length of the sequence. The 
base case is trivial as /o.-i = {0,0} and fcfl are dummy, 
c S [0, z^]. Suppose the statement holds for j — 1 variables. 

Consider the variable Xj. Suppose, by contradiction, that 



Ui < /< 



Qc+l,j-! I-Cj 



Ir.., < I 



Then either qcj < qc+ij or qcj 



c+i-r 



By induction hypothesis, we know 



that /cj_i > /c+ij-i, hence, either qcj-i > Qc+ij-i or 

1c,j-l — <lc+l,j-l,lc,j-l > ^c+lj-1- 

We consider three cases depending on whether Xj is a pe- 
nalizing variable, an undetermined variable or a neutral vari- 
able. 

Case 1. Consider the case Xj G Pk- If qc.j-i = oo then 
Qc+i.j-i = oo by the induction hypothesis. Hence, by Algo- 
rithmls] linel9] fcj and fc+i.j are dummy and equal. Suppose 
Qcj-i 7^ OO. Then we consider four cases based on relative 
values of qc,j' , qc+i,j' ,lc,r , lc+i,r , f e {j - l,j}. 



Case la. Suppose qc.j < qc+i,j and qc,j-i > qc+i.j-i 



By Algorithm |3] lines |9] and [TO 

9c+ij < qc 

Qc+lJ-l < 
diction. 






l.nHence 
Qc+ij 1 



9c J > Qc,j-i and 
Qc,j < qc+i.j imphes 
1. We derive a contra- 



• Case lb. 



and 



ge+l ;- li hj-l ^ 'c+lj' 



10 



and 
Hence 

9cj < 



qc.j 

9cj 

Qc+l.j 



< Qc+ij ana g^j- _ 
1. By Algorithm [3] lines |9] 

< 



qcj then lf,,j 



Suppose qcj 
> I 
> qc,j-i and qc+i,j < Qc+Lj-i -f 
< 9c+ij implies Qc+ij-i = Qc^j-i - 
< qc+ij-i + 1- Hence, qc+ij-i = 

len) by Algorithm [3] Une |9] As 

1 then Ic+ij-i e {len, oo} by 



9cj-i = Qcj and qc+ij = Qc+ij-i + l-_As qc._ 



1 e [1 

1c+i,j — Qc+ij-i 4 

Algorithm [3] line [TO] This leads to a contradiction as 

c+lj-l- 

Suppose gcj = qc+i,j,lc3 < 



Case Ic. 



Qc.j-i > Qc+ij-i- Symmetric to Case lb. 
Case Id. Suppose g^j — qc+ijJc.j < 



lc+i,j and 



> I 



tc+lj 



and 



c+lj-l 



By Algorithm 



-iJc.j-l 

qc,] > 9cj-i and qc+i,j < Qc+ij-i 4 
: qc+i,j imphes Qc+ij-i = gcj-i 



< 



< 



Qc+i j-i + 1. Therefore, either qcj 



-1 or qcj = q, 



cj-l 



+ 1 A 



Qcj-l — Q c+l j 
lines[9]and[TO 
Hence, qcj - 
Qcj = Qc+ij 

IcJ-l A Qc+lJ = Qc+lJ 

Qc+ij — Qc+ij-i + 1- 

If qc,j = Qcj-i and ^c+ij = Qc+ij-i then Lj^i G 
[1, len) and Ic+ij-i G [1, Zen) by Algorithm [3]line [9] 
Hence, Icj — Icj-i + 1 and Ic+ij = Ic+ij-i + 1. As 
Icj-i > lc+ij~-i, then Icj > I'c+ij- This leads to a 
contradiction with the assumption Icj < Ic+ij- 
If ^cj = gcj-i + 1 A qc+ij = qc+i,j-i + 1 then 
Icj-i G {len, oo} and Ic+ij-i G {len, oo} by Al- 
gorithm [5] line [TO] Hence, Icj = 1 and Ic+ij — 1- 
This leads to a contradiction with the assumption Icj < 

Ic+lJ- 

Case 2. Consider the case Xj G Uk- If Qcj-i = oo then 
^i.j_i = cx) by the induction hypothesis. Hence, by Algo- 

rithm[3] line 12 fcj and fc+ij are dummy and equal. 

Suppose qcj-i 7^ oo. Then we consider four cases based 

on relative values of qcj',qc+i,j',lc.3',lc+ij', f e {j - 

Case 2a Suppose qcj < qc±ij and qcj-i > 

we know 
9cj-i < 



q, 



qc+ij-i- By Algorithm|3] lines [T2] and 
that qc+ij-i < Qc+ij < qc.j-i an( 



< 



qc-ij-i' 

< 

< 



By induction hypothesis, Qc+ij- 
Hence, if qcj 

qc+ij < Qcj-i- 



< 



9cj 

9cj-i < qc-ij-i- Hence, if qc,j < Qc+ij then 
Qcj-i < gcj < qc+ij < Qcj-i- Therefore, if 
qcj < qc+ij then we derive a contradiction. 

Case 2b. Identical to Case 2b. 



Suppose qcj = 

<Zc+ij"-i- As 

= qc+1 



> 



Ic+ij and 



-I J ( line[T2]i 
< qcj-i Wr 



qc+ij 7 i-cj 

9cj-i 7^ 9c+ij-i then 
We also know qcj-i < 



Case 2c. 

9c J -1 > 
9c+ij-i = 

9c,j < qc+ij < Qcj-i Irom Case la. Putting every- 
thing together, we get Qcj-i < ^cj < qc+ij-i < 
qcj-i- This leads to a contradiction. 

Case 2d. Suppose qcj = qc+ij,lc 



qcj-i — qc+ij-i,lcj-i 



>l 



'cj 



< I 



c+l,J 



and 



c+l,i- 



_ 1 . As we know from 



Case la Qc+ij-i < Qc+i^ < Qcj-i, Qcj-i < Qcj < 
9c-ij-i and Qcj-i < Qcj < Qc+ij < Qcj-i- Hence, 
Qc+i,j-i — <lc+i,j = Qcj-i ~ Qc.j- 
Consider two subcases. Suppose Qcj-i < Qc-ij-i- 
Then L 



icj — oo (line 12 1. Hence, our assumption 
Ir-i < lc-t-T--i is false. 



Suppose Qcj-i 



gc-ij-i- 



lc,j — oo (line 12 1. Hence, our assumption Icj < I 



If ^f^-ij-i = len then 



'c+lj" 



is false. Therefore, /c-i.j-i G [1,/en) and lc,j-i — 
Ic-i.j-i + 1- By induction hypothesis as Qc+ij-i = 
Qcj-i = Qc-Lj-i then /c+ij-i < ^cj-i < ^c-ij-i- 
Hence, ^cj-i G [li^c-ij-i] ^ [1,/en). Therefore, 
lc+i,j — ^c,j-i + 1 < 'c-ij-1 + 1 = hj-i- This 
contradicts our assumption Zcj < ^c+i.j- 
Case 3. Consider the case Xj e A^^. This case follows 
immediately from Algorithm |3] line [14] and the induc- 
tion hypothesis. 

n 



F Proof of Lemma |6] 

Proof of correctness. We prove by induction on the length of 
the sequence. Given fcj we can reconstruct a corresponding 
set of sequences Scj by traversing the table backward. 

The base case is trivial as xi G Pk, /o,o — {l; 1} and 
/c,o = {oo, oo}. Suppose the statement holds for j — 1 vari- 
ables. 

Case 1. Consider the case Xj G Pk- Note, that the cost can 
not be increased on seeing Xj e P^ as cost only depends on 
covered undetermined variables. By the induction hypothe- 
sis, S'cj-i satisfies conditions [T]-p] The only way to obtain 



Scj from 5c' j-i, c' € [0, z^, ], is to extend last{Sc, 



i-i 



to 



cover Xj or start a new sequence if \last{Sc,j-i)\ = len. If 
5c J- 1 does not exist then 5c. idoes not exist. The algorithm 
performs this extension (lines|9|and[T0|. Hence, Scj satisfies 
conditions [T]-2] 

Case 2. Consider the case Xj £ Uk. In this case, there exist 
two options to obtain Sc.j from from 5c' j-i, c' S [0, z^]. 

The first option is to cover Xj . Hence, we need to extend 
last{Sc-i.j~i). Note that we should not start a new sequence 
if last{Sc-ij^i) = len as it is never optimal to start a se- 
quence on seeing a neutral variable. 

The second option is not to cover Xj. Hence, we need to 
interrupt last{Scj-i). 

By Lemma [5] we know that fc.j-i < fc-i.j-i, < c < 
C By the induction hypothesis, 5c.j-i and 5c- ij-i satisfy 
conditions [Tl-H] Hence, 5c.j-i < 5c- i.j-i. 

Consider two cases. Suppose 1 5c j-i| < |5c-ij-i|. In 
this case, it is optimal to interrupt last{Scj^i). 

Suppose |5c,j-i| = |5c-ij-i| and \last{Sc,j^i)\ < 
\last{Sc-i.j^i)\. If \last{Sc-ij-i)\ < Zen then it is optimal 
to extend last{Sc-i,j-i). If \last{Sc-i,j-i)\ — len then it 
is optimal to interrupt last{Sc.j-i), otherwise we would have 
to start a new sequence to cover an undetermined variable Xj, 
which is never optimal. If 5c.j-i and 5c-ij-i do not exist 
then Scj does not exist. If 5c.j-i does not exist then case 
analysis is similar to the analysis above. 



This case-based analysis is exactly what Algorithm [3] does 
in lines 12 and 12 Hence, Scj satisfies conditions [TlBr 



Case 3. Consider the case Xj € Nk. Note that the cost can 
not be increased on seeing Xj E N^ as cost only depends on 
covered undetermined variables. By the induction hypothe- 



sis, 5c., -1 satisfies conditions 



Sc,j from 5, 



c'J-l 



d e [0, 



4] The only way to obtain 
is to interrupt last{Sc,j-i). 



If 5c J- 1 does not exist then Sc.j does not exist. The algo- 
rithm performs this extension in line 14 Hence, 5c j satisfies 
conditions [TJHl D 

Proof of complexity. The time complexity of the algorithm is 
0(nmax(zc)) = 0{n^) as we have 0(nmax(2:c)) elements 
in the table and we only need to inspect a constant number of 
elements to compute /(c,j). D 

G Proof of Lemma H] 

First, we present explicitly the algorithm for detecting disen- 
tailment. 



Algorithm 4: WEiGHTEDSPRiNGYFocus(a;o, . . . ,x„-i) 

1 for c G — 1..2:/' do 

2 for j e — l..n — 1 do 

3 fc.j ^ {00,00,00}; 

4 /o,-i ^ {0,0,0}; 

5 for j e 0..n — 1 do 
f. for c e 0..j do 

7 \ixj G Pfc then /* penalizing */ 

8 if (^c,j-l ^ [1, len)) V (gc,j-l — oo) then 

9 id <— {9c, 3-1, Ic, j-l + 1, 'ic.j-l}; 

10 else 

11 /c,j ^ {9c,j-i + 1,1,0}; 

12 \^Xj G Uk then /* undetermined */ 

13 if (i<:-l,3-l S [1, ien) A (jc-l,j-l = 
<Jc,3-i) V ((3c,j-i = 00) then 

W /c,J •*- {gc-l,j-l, la-l,j-l + 1, /lc-1,3-1} 

15 else 

16 /c,j <- {(Jc,3-1,00, 00} 

17 if a;^ tE N^ then /* neutral */ 

18 if (/c-l,j-l ^ [1; ^en) A /ic-i.j-i ^ 

[l,/i) A gc-1,3-1 = qa, j-i) V (gc,j-i = 00) then 

19 /c,j <- {?c-l,j-l, la-l,j-l + 1, /ic-l.i-l + 1} 

20 else 

21 /c.j •*- {?c,3-l,00, 00} 



22 return /; 



The algorithm is based on a dynamic program. Each 
cell in the dynamic programming table fc.j, c G [0, Zc'^], 
j e {0, 1, ... ,n — 1}, where z^ = max{zc) — \Pk\, is a 
triple of values gc.j, Icj and hcj, fcj = {qcj,lcj,hc,j}, 
stores information about Scj . The new parameter hcj stores 
the number of neutral variables covered by last{Sc,j). Algo- 
rithmfflshows pseudocode of the algorithm. 

The intuition behind the algorithm is as follows. We em- 
phasize again that by cost of a cover we mean the number of 
covered undetermined and neutral variables. 

• If Xj G Pk then we do not increase the cost of Scj com- 
pared to 5c,j-i as the cost only depends on Xj G Uk- 
Hence, the best move for us is to extend last{Sc,j-i) or 
we start a new sequence if it is possible. This is encoded 
in hnes[9]and 11 of the algorithm. 



If Xj G Uk then we have two options. We can ob- 
tain Scj from Sc-i.j-i by increasing cst{Sc-ij-i) by 
one. This means that Xi will be covered by last{Scj). 
Note that this does not increase the number of covered 
neutral variables by last{Scj) as we can always set 
Xj = V, V > k. Alternatively, from Scj-i by interrupt- 
ing last{Sc,j-i) if necessary. This is encoded in lines 14 
and [16] of the algorithm. 

If Xj e Nk then we have two options. We can obtain 
Scj from Sc-ij-i by increasing cst{Sc-i.j-i) by one 
and increasing the number of covered neutral variables 
by last{Sc,j-i)- Alternatively, from Scj-^i by inter- 



rupting last{Scj-i) (lines 19-21 1. 

The proof of correctness mimics the corresponding proof 
for the WeightedFocus constraint. We can now prove 
Lemma [8] 

Proof. The main idea is identical to the proof of the 
WeightedFocus constraint. We only highhght the dif- 
ferences between the WeightedFocus constraint and the 
WeightedSpringyFocus constraint. 

Consider a variable-value pair Xi — v, v > k. The 
only difference is in the fourth option. We denote h{sij) 
the number of neutral variables covered by Sij. Similarly, 

• The fourth and the cheapest option is to glue 

last{Sci,i-i), Xy and last(Sc2,n-i-2) to a single se- 
quence if \last{Sci.i-i)\ + \last{Sc2.n-i-2)\ < ^enand 
h{last{Sci,i-i)) + h{last{Sc2,n-i-2)) < h. Hence, 

'-'ci.z-l ~ 'Jci.i-l \ tasi(Dci,i_l), '-'c2,n-'i-2 = 

Sc2,n-i-2 \ lO'St{Sc2,n-i~-2) and s' is a concatenation 
of last{Sci,i~i),x = V and last{Sc2,n-i~2)]- Then the 
union S — S'^^ ^_i U S'^^_^_^_2 ^ {^'} forms a cover: 

CSt{S) = C1+C2 + I, |5| = |S'ci,i-l| + |5'c2,n-i-2| -1 

and h{S) = h{last{Sci,i-i)) + h{last{Sc2.n-i-2))- 

The rest of the proof is analogous to WeightedFocus. 

Consider a variable-value pair Xi = v, v < k. The main 
difference is that we have the second option to build a support. 
Namely, we glue SciA^i, Xi and Sc2,n~i-2- Hence, if ci + 

C2 + 1 < z^ , \last{Sc^^i_i)\ + \last{Sc2,n-i-2)\ < len 
and h{last{Sci.i-i)) + h{last{Sc2,n-i-2)) < h then we can 
build a support for Xi = v. The rest of the proof is analogous 
to WeightedFocus. D 



